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The computation of the critical exponent rj characterizing the universal elastic behavior of crys¬ 
talline membranes in the flat phase continues to represent challenges to theorists as well as computer 
simulators that manifest themselves in a considerable spread of numerical results for rj published 
in the literature. We present new insight to this problem that results from combining Wilson’s 
momentum shell renormalization group method with the power of modern computer simulations 
based on the Fourier Monte Carlo algorithm. After discussing the ideas and difficulties underlying 
this combined scheme, we present a calculation of the renormalization group flow of the effective 2d 
Young modulus for momentum shells of different thickness. Extrapolation to infinite shell thickness 
allows to produce results in reasonable agreement with those obtained by functional renormaliza¬ 
tion group or by Fourier Monte Carlo simulations in combination with finite size scaling. Moreover, 
our new method allows for the first time to obtain a decent estimate for the value of the Wegner 
exponent oj that determines the leading correction to scaling, which in turn allows to refine our 
numerical estimate for rj previously obtained from precise finite size scaling data. 
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INTRODUCTION 

The paradigm of the renormalization group (RG) 
is without doubt a cornerstone of modern theoretical 
physics with countless applications, and has been enor¬ 
mously influential in many areas of science beyond its 
origins rooted in high energy physics and statistical me¬ 
chanics [1]. Indeed, the abstract RG concept may be 
regarded as a cleverly organized successive divide-and- 
conquer strategy to deal with problems that involve a 
large number of mutually coupled degrees of freedom. 
Yet, concrete applications of an RG scheme may superfi¬ 
cially appear to look very different from one another. In 
the present article we shall concentrate on Wilson’s mo¬ 
mentum shell RG (MSRG) approach to the field-theoretic 
formulation of critical phenomena at second order phase 
transitions [2] . The MSRG is certainly an invaluable con¬ 
ceptual tool both for abstract reasoning as well as in a 
first qualitative or even semi-quantitative analysis of a 
given problem. In a nutshell, one writes the underly¬ 
ing Hamiltonian in terms of Fourier amplitudes /(fc) of 
the underlying helds. Imposing a wave vector cutoff A, 
one tries to identify an effective Hamiltonian as it would 
emerge after having integrated out all microscopic de¬ 
grees of freedom that describe the physics of the system 
below scales of size 1/A. In this effective Hamiltonian 
only those couplings are kept that are regarded as impor¬ 
tant in the long wavelength limit, while the effect of all 
other couplings that are related to the eliminated short- 
ranged degrees of freedom is absorbed into an assumed 
renormalization of these surviving couplings. The fact 
that the choice of the cutoff A is arbitrary suggest to 
iterate this prescription as follows. The effects of the 
“fastest” degrees of freedom which reside in a momen¬ 
tum shell A/6<|fc|<A, 5>1 beneath the cutoff A are 
successively integrated out from the partition function. 


which gives rise to a yet another set of modified coupling 
constants. On properly rescaling lengths and “renormal¬ 
izing” the held, one derives a how pattern in the space of 
coupling constants. An analysis of the hxed points (FPs) 
of this how then allows to explain the phenomenon of 
universality and to extract numerical values for the criti¬ 
cal exponents. Unfortunately, concrete analytical imple¬ 
mentations of this program usually rely on some type of 
perturbative approximation, and calculations frequently 
become intractable beyond one loop order. Thus, for 
actual numerical calculations other approaches like the 
held-theoretic RG [3] or the functional RG [4] are pre¬ 
ferred, or one resorts to real space computer simulations 
in combination with hnite size scaling (FSS) 0(0]. 

For the task of implementing Wilson’s MSRG scheme 
in a simulation, real space MG approaches are obviously 
not very well suited. On the other hand, our Fourier 
Monte Garlo algorithm (FMG) [ZHII] is tailor-made for 
this problem. Recently we have demonstrated that it 
allows to follow the MSRG prescription step by step in 
simulation [alia. This is quite appealing, as it elim¬ 
inates the perturbative approximations and the under¬ 
lying need for a “small parameter” from the concrete 
application of the MSRG, thus representing a truly non- 
perturbative implementation. On the other hand, like 
in any humanly possible MSRG calculation, one is still 
forced to project the calculated RG flow from its na¬ 
tive infinite-dimensional coupling space to a suitable low¬ 
dimensional subspace spanned by a finite number of effec¬ 
tive coupling parameters. Even though no perturbative 
approximation is involved, the presence of this inevitable 
projection, which amounts to ignoring the effects of the 
remaining directions in the space of coupling constants, 
may do substantial harm to the achieved numerical pre¬ 
cision. 

Inspired by early analytic work by Bruce, Droz and 


2 


Aharony [H, subsequent work [niin] indicates that by 
optimizing the results with respect to the parameter b 
that governs the thickness of the momentum shell (we will 
discuss below in more detail how this works) MSRG can 
indeed be turned from a qualitatively to a quantitatively 
useful tool. Up to date this has only been demonstrated 
for a particularly convenient model system, namely the 
long-ranged Ising model of Fisher, Ma and Nickel [H]. 
One purpose of the present paper is to test the ideas 
put forward in Refs. mun] on a nontrivial real world 
problem, namely the elastic behavior of crystalline mem¬ 
branes in the flat phase. It is well known m that the 
corresponding universal behavior of long wavelength fluc¬ 
tuations is governed by single exponent rj. However, as 
noted in Ref. [T5], a glance at the existing literature re¬ 
veals a considerable spread of numerical results for 77, 
obtained from a variety of analytical approaches like e.g. 
self-consistent field theory [T9l|20], e-expansion [21], large 
d expansion |22| and functional RG [23l[2^ , or from sim¬ 
ulation approaches derived in real space (see e.g. Refs. 
[ISIIIT] ). In our own simulations (Ref. [IH]), which are 
based on our Fourier MG algorithm in combination with 
FSS, we have noticed a rather strong influence of cor¬ 
rections to scaling, which indicates the importance of 
properly taking into account the role of RG-irrelevant 
couplings if one aims at high numerical precision. In¬ 
terestingly, up to date nobody seems to have succeeded 
in deriving a numerical estimate of the exponent uj gov¬ 
erning the corrections to scaling (cf. [28|) of a crystalline 
membrane in its flat phase. It is the second goal of the 
present paper to provide such a numerical estimate. 


A DIFFERENT VIEW ON STANDARD MSRG 

As explained in the introduction, MSRG is a fairly 
standard method. Yet, for the convenience of the reader 
we will summarize the main steps, emphasizing those as¬ 
pects that are particularly relevant to our present ap¬ 
proach. A MSRG transformation can be performed on 
an arbitrary Hamiltonian formulated in terms of 

the Fourier amplitudes f{k) of a field f(x) defined for 
wavevectors of moduli up to a chosen momentum space 
cutoff A. The formal vector K holds all “coupling con¬ 
stants” that are admissible for the symmetry constraints 
imposed on the underlying system. Let K. denote the 
infinite-dimensional space of all such coupling vectors. 
One chooses a shell thickness parameter 6 > 1 and splits 
the Fourier amplitudes /(fc) = /<(fe) + />(fe) into “slow” 
and “fast” contributions 

/<(fc) = 0(A/6-|fe|)/>), (1) 

/>(fc) = 0(|fc|-A/6)/>), (2) 


where 0 denotes the Heaviside step function. Functional 
integration over the fast modes 

(3) 

then yields a new coarse-grained Hamiltonian for remain¬ 
ing slow modes and induces a mapping K —>■ K. This 
coarse-graining step is followed by a rescaling k' = bk of 
“momenta” (i.e. inverse length) scales and restoration of 
the original cutoff A. Finally, one performs a renormal¬ 
ization 

h{k'/b) = z{b,K)f{k') (4) 

of field amplitudes with 

z{b,K)=b^-^f'^-'^ . (5) 

Here d is the spatial dimension, [/] is the canonical mo¬ 
mentum dimension of /, and the so-called anomalous 
dimension ri[K]/2 is a function which characterizes the 
specific FP to be investigated (see below). For the ex¬ 
ample of a coupling constant C multiplying a monomial 
containing n powers of the field / and p spatial deriva¬ 
tives in the effective Hamiltonian, it is straightforward 
to show that the coarse-grained coefficient C undergoes 
a total rescaling 

C" = (6) 

where p is implicitly accounted for in the canonical mo¬ 
mentum dimension [C] = d — n ■ [f] — p. 

Consecutive application of these three steps induces a 
mapping K —>■ K', which defines the RG transformation 
TZb : /C —>■ /C. The crux of the whole construction is the 
observation that as a result of the rescaling operation, the 
correlation lengths of systems at K and K' are related by 
£,[K'] = ^[K]/b. At a FP RT* = UbiK*) this leaves only 
the possibility of an infinite or zero correlation length. 
Each such FP characterizes a different universality class 
of critical behavior, and nontrivial behavior is, of course, 
found for infinite correlation length. 

In principle, the operation TZb can be defined for any 
6 > 1, and satisfies the eponymous semi-group property 

Ti-bib2 = T^bi °T^b2, (V 

which is paramount to the emergence of power laws that 
dominate the subsequent analysis as well as to the inde¬ 
pendence of the associated exponent values of the par¬ 
ticular choice of the shell thickness parameter b. In the 
vicinity of K* where TZb can be linearized, most direc¬ 
tions in the space K, turn out to be exponentially at¬ 
tractive (“irrelevant”), while typically only one or two 
are exponentially repulsive (“relevant”), and thus must 
be carefully tuned to “reach” the FP K* under succes¬ 
sive iteration of TZb- Ultimately, this explains the ob¬ 
served universality of critical phenomena. The RG flow 
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resulting from the above scheme is defined in the infinite¬ 
dimensional coupling constant space JC. In practical cal¬ 
culations, one is nevertheless forced to limit ourselves to 
working with effective Hamiltonians, i.e. Hamiltonians 
parametrized exclusively by coupling 
vectors KeS € tCeS taken from a low-dimensional linear 
subspace ICeS C /C of dimension, say, doff) spanned by 
the relevant and the least irrelevant directions w.r.t. the 
FP K*. In terms of suitably chosen coordinates in the 
space 1C, the projection : /C —)■ /Ceff onto this finite¬ 
dimensional space assumes the form 

7reff(iFi,iF2,..-) = (8) 

Except for trivial cases, ICeS is not an invariant subspace 
under the action of TZb i.e. TZb and TTeff do not com- 
mnte, because “new” couplings are inevitably generated 
from a generic effective Hamiltonian under the coarse- 
graining operation, regardless of our ability to perform 
the coarse graining operation exactly or by some approx¬ 
imate method. If the subspace /Ceff C /C has been chosen 
properly, the “missing” directions will only correspond 
to strongly irrelevant directions in coupling space, whose 
influence will be exponentially suppressed. In mathe¬ 
matical terms, any humanly possible MSRG calculation 
amounts to replacing the exact RG transformation TZ by 
the effective transformation 

TV^b ■= TTeff oTlbO TTeff ■ (9) 

The crucial observation is, however, that, in contrast 
to TZb the effective transformations do not strictly 
form a half-group, i.e. 

TZ^^b.b, + o TV^b^ , (10) 

since in the composite operation on the right side the 
additional irrelevant couplings generated by 7?.®®;, will be 
“lost” in the subsequent application of TTeff. Of course, 
nothing can prevent us from studying iterations of the 
map TZ^^b in a manner similar to TZb- The harm that the 
failure of to close under composition causes to the 
subsequent analysis depends on the “production rate” of 
coupling components generated during the coarse grain¬ 
ing step that fall outside of /Ceff, which in turn is con¬ 
trolled by the shell thickness parameter b. In particular, 
the projection TTeff (IT*) of the “true” infinite-dimensional 
FP K* of TZb does generally not produce a FP of 7T®®;,. 
Instead, the locations of FPs iTes* = lTeff*(6) € /Ces of 
the transformations 7C®®b will generally be b-dependent. 
Furthermore, the same is true for the numerical values of 
critical exponents calculated from a linearization of 7T®®b 
around Kee*{b). 

In summary, even though the results of the exact RG 
prescription in infinite-dimensional coupling space /C are 
guaranteed to be independent of the arbitrary parameter 
b, the projection TTeff to the low-dimensional space /Ceff 
introduces such a ^-dependence that encodes the effects 


of the remaining irrelevant directions. While this seems 
to look pathological at first sight, it actually allows to 
optimize the resulting calculation scheme by determin¬ 
ing the value b* at which the drift of Kgff{b) becomes 
stationary. In this respect, our philosophy is similar to 
that of other approaches in which an arbitrary parameter 
is introduced whose value would drop out of the results 
of exact theory but nevertheless may be used to optimize 
an approximated version. A nice example illustrating the 
power of such a strategy is H. Kleinert’s “variational per¬ 
turbation theory” |29j . However, it is very important to 
keep in mind that the present 6-related “pathologies” are 
non-perturbative in the sense that they do not originate 
from the use of any perturbative approximation in evalu¬ 
ating the CG step, but purely arise from the necessity to 
limit ourselves to considering a finite number of couplings 
in a real world calculation. Amusingly, these effects are 
neither noticed in standard perturbative MSRG calcula¬ 
tions, where it is extremely convenient to consider mo¬ 
mentum shells that are infinitesimally thin, since in the 
limit Ab := 6 — 1 —>■ 0+ the appearing Feynman integrals 
are usually much easier to evaluate than for finite Ab, nor 
in most popular real-space RG schemes where the value 
of b is usually dictated by the decimation scheme cho¬ 
sen for the given lattice topology. In fact, it is difficult to 
find any papers that use momentum shells of finite thick¬ 
ness for anything beyond qualitative arguments. One no¬ 
table exception is the work of Bruce, Droz and Aharony 
[14] , who argued that the influence of irrelevant couplings 
in perturbative calculations of the exponents of a stan¬ 
dard short-ranged Landau-Ginzburg (LG) model should 
be greatly diminished in the limit of large b. And indeed, 
notice that 6* —>■ oo and 6* —)■ 1 are the only values of b* 
that allow to reconcile the expected 6-dependent features 
discussed above with the validity of the usual semi-group 
property 7Tf,.2 = TZb* ° TZb* ■ 

Our recently developed FMC method is non- 
perturbative by definition and necessarily uses momen¬ 
tum shells of finite thickness, since our simulations 
are done for a finite lattice of linear size L with lat¬ 
tice constant a = 1, which implies a minimum spac¬ 
ing of Aki = 2Tr/L between components of adjacent 
wavevectors. Thus, it is perfectly suited to study the b- 
dependence of TZ^^b and check the predictions of Bruce et 
al. that had been derived with the use of the e-expansion. 
Of course, due to the discrete nature of the Brillouin 
zones of our finite systems neither the limit 6 —> 1 nor 
the limit 6 —>■ oo are directly accessible, but we can mon¬ 
itor or even try to extrapolate the behavior of the corre¬ 
sponding observables towards these limits. 

For the purpose of putting our ideas to the test, the 
short range LG model used in Ref. [14] is not very suit¬ 
able in view of the numerical smallness of its exponent 
T] = e^/54-|-0(e^). Instead, in Ref. [T3| we considered the 
long-range generalization of the LG model introduced by 
Fisher, Ma and Nickel in Ref. m- This model was par- 
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ticularly convenient since the exponent r] of its Wilson- 
Fisher FP is exactly known, thus saving the numerical 
effort to determine it numerically from the simulation 
data. In addition, detailed analytical calculations and 
quite precise Monte Carlo data were available for compar¬ 
ison [3D]. Using our FMC implementation of MSRG, we 
were indeed able to observe the ^-dependence of KeB*{b) 
and its associated exponents. However, contrary to our 
initial expectations, it turned out that the best accuracy 
was not obtained in the large b limit. Instead, for varying 
b the FP b i— Kes*{b) moves along a “trajectory” in the 
plane /Ceff that exhibits a turning point at a certain shell 
thickness b* that was actually found to be rather close 
but distinct from 6 = 1, and for this distinguished value 
b* we observed that the values of the critical exponents 
v and Lo were in excellent agreement with the benchmark 
results derived in Ref. |3D|. A systematic study for dif¬ 
ferent system sizes revealed the surprising discovery that 
1 < 6* < oo is not a finite size effect. Nevertheless, we 
speculate that this peculiar finding is highly specific to 
the model of Fisher, Ma and Nickel, and we still expect 
that usually 6* —)■ 1 or 6* —>■ oo will instead be found 
in other systems. The rest of the paper will therefore 
be devoted to the application of our ideas to a real-world 
system, whose critical properties are still an active area of 
research: the elastic behavior of crystalline membranes. 


FMC IMPLEMENTATION OF MSRG FOR 
CRYSTALLINE MEMBRANES 

As explained in detail in Refs. [17l|3T], the flat phase of 
a crystalline membrane is conveniently described in the 
so-called Monge parametrization, which amounts to spec¬ 
ifying a scalar “height” function f{x) that measures the 
out-of-plane deformations of the membrane with respect 
to a two-dimensional reference plane, which we take to 
be of size L x L with periodic boundary conditions un¬ 
derstood. The long-wavelength physics of the system is 
captured by the Fourier modes 

f{q) = 0iA - |q|) J , (II) 


Its first contribution, the bending energy, is represented 
by a local dispersion term as found in a standard LG 
model, except that the usual gradient term (V/)^ is re¬ 
placed by a Laplacian (A/)^. In addition, while its sec¬ 
ond contribution also involves four powers of /, the non¬ 
local characteristic of the generalized convolution (12) 
hints at physics that is quite different from that of the 
standard LG model. As discussed in Refs. [HI Eg, this 
non-locality encodes an effective long-range anharmonic 
self-interaction of the out-of-plane deformations / me¬ 
diated by in-plane phonons that had been integrated 
out in the calculation steps leading to (13). H)/] in¬ 

volves only two coupling constants, namely the bend¬ 
ing stiffness k and the the effective 2d Young modulus 
K = 4/r(/r -|- A)/(2^ -|- A) composed from the in-plane 
Lame constants A and /r of the membrane. Implicit in all 
these constants as well as in the formulas (13 but sup¬ 
pressed in our present notation is a dependence on the 
cutoff A. 

To implement our FMG algorithm, we replace the 
membrane reference plane by a square L x L lattice with 
N = LF' sites and lattice constant a = 1, and keep the 
imposed periodic boundary conditions. Assuming with¬ 
out loss of generality L to be even, we may parametrize 
wavevectors inside the full first Brillouin zone of this lat¬ 


tice by Qi = 27rTOi/L, to, = —L/2-I-1,. .., 0, ..., L/2, and 
the above integrals over the Brillouin zone are replaced 
by finite sums. In view of the rectangular structure of 
the underlying lattice, it is natural to replace spherical 
cutoffs A that are convenient in analytic continuum cal¬ 
culations by a more suitable cubic version. Parametrized 
by an integer I, in our simulation a cutoff A = 2-^1 /L is 
applied to each separate wavevector component, and in 
order to avoid problems with “umklapp” terms and min¬ 
imize effects of lattice anisotropy, it is recommended to 
choose I L/ 2. To implement the coarse graining step in 
FMG, we furthermore choose an inner cutoff A' = 2^1' /L 
with 0 < I' < 1. The shell thickness parameter is then 

given by 6 = A/A' = 

The discrete Fourier transform convention 




where the Heaviside step function is used to impose a 
cutoff A in the space of wavevectors. Formally embedding 
the vectors q,Q in and abbreviating Q = Q/\Q\, we 
define 

HQ) = / ^ (Q X 9 ) ' - 9 ) • (12) 

In terms of this generalized convolution, the effective 
Hamiltonian that describes the universal properties of 
the flat phase at long wavelengths is the given by 


with inversion 




(15) 


l9d<A 


in which the Fourier amplitudes are extensive quantities, 
may look somewhat asymmetric, it proves to be conve¬ 
nient in comparing discrete to continuous formulas. Since 
the membrane’s elastic free energy does not depend on 
the average distance of the membrane to the Monge ref¬ 
erence plane but merely on variations of its height, only 
derivatives of / enter in the in the continuum formula¬ 


tion (13). Therefore we can further assume without loss 


of generality that /(O) = 0. In terms of these discrete 
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amplitudes, the above formulas (12) and (13) are replaced 
by 

•^(Q)=E(Q''9)"/(<7)/(Q-q) (16) 


and 


n^if] = ^ E ^ E ’ (17) 

q^O Q^O 


where 


K K 

= Kn = ^. 


(18) 


In view of the extensive discussions already available 
in the literature (cf. Refs. 151 fTOl 1551 and [33]) and the de¬ 
tailed layout of the specihc implementation for the case 
of crystalline membranes presented in the companion pa¬ 
per [T8| , we would like to keep the description of the basic 
Fourier Monte Carlo algorithm and its general properties 
at a minimum in the present paper. However, it turns 
out that setting up the coarse graining step of MSRG for 


Q 


with 

-F(Q) = ^ (q X q) ' f{q)f{Q - q) , (21) 

q 


from we wish to extract the two CG relations Kjq i—>■ kjq 
and Km ^ Km, i.e. k k and K ^ K. For this 
purpose, we determine the value of the CG Hamiltonian 


(20) by restricting the MG sampling to certain “tracer 


configurations” defined by a particularly simple and con¬ 
venient choice of their slow mode parts. In terms of sim¬ 
plicity, our preferred type of such a tracer configuration 
would certainly be that of an isolated “dumbbell” of just 
two slow modes with a common uniform real-valued am¬ 
plitude at the fixed wave vector ±fc. This dumbbell is 
surrounded by the shell of nonzero fast modes, but all 
other slow modes are put to zero. In formal terms, the 
slow parts of such dumbbell tracer configurations (q) 
dehned with respect to ±fc are restricted to be of type 


/4'')(q) = /d(Vfc+W). 


( 22 ) 


As explained in detail in Refs. [3 m, for this class of 
tracer configurations one now performs a multicanoni- 
cal type of simulation of e.g. the Wang-Landau type, in 
which the probability distribution P{fd) of the ’’reac¬ 
tion coordinate” fd in the “bath” of fast modes is cal¬ 
culated. A polynomial fit of —\nP{fd) then yields a 


a crystalline membrane requires to define different MG 
moves for slow and fast modes of the so-called “tracer” 
configurations to be dehned below. In the standard cu¬ 
bic FMG scheme the momentum shell corresponding to 
a prescribed pair of cutoffs A' < A is, of course, dehned 
as the set of wave vectors with components pi,i = 1,2 
subject to the constraints |pi| < A and max^ \pi\ > A'. 
MG move of fast modes f>{q) are performed by picking a 
random wave vector p from this shell, choosing a random 
complex number inside a circle |e| < p of radius p around 
0 in the complex plane, and considering the shift 


f{.q)^f{q)p^5q^p + (*5q-p. (19) 


Taking advantage of the special convoluted structure (12) 
of the anharmonic term appearing in (13), it is then pos¬ 
sible to calculate the resulting change in energy in an 
efficient way, as is explained in detail in Ref. [18]. 


Integrating out these fast modes by means of an FMG 
simulation should then produce a coarse-grained Hamil¬ 
tonian of general structure 


E 

Q 


\P{Q)? + 0{f) 


( 20 ) 


set harmonic and lowest order anharmonic coefficients 
a 2 {k), a 4 (fc),... for each chosen wave vector k. Gompar- 
ison of these coefficients with the general fc-dependent 
structure of the bare effective Hamiltonian then allows 
to determine a “new” set of bare parameters. In other 
words, one obtains all the information required for com¬ 
pleting the coarse graining step of the MSRG prescrip¬ 
tion. 

For LG type of models with a local anharmonic en¬ 
ergy contribution, this class of tracer configurations al¬ 
lows to determine the flow of coupling parameters. Un¬ 
fortunately, however, for our present problem the dumb¬ 


bell class (22) is insufficient to capture the flow of the 


anharmonic part of the bare Hamiltonian. In fact, in the 
formula 


^^<\Q) = (Q X kf [5Q.2k + 25q + 5Q+2k\ fd (23) 


that results for amplitudes P{Q) built exclusively from 
the slow part /<^^(q) = fd{5q~k + 5q+k), the vector Q 
is constrained to be either zero (which is forbidden) or 
parallel to fc, in which case the leading cross product 
vanishes, i.e. P^^\Q) = 0. To overcome this difficulty, 
we instead consider “cross” tracer configuration with slow 
parts of type 

/<(?) = /c fq^k + 5q-k^ + 5qj-k + 5q+kf , (24) 
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where /c € K and |fc| = fe • fe'*' = 0, which map out 
a symmetric “cross” spanned by two orthogonal vectors 
k and fc"*- of equal length around 0 with one common 


real-valued amplitude, all remaining slow modes being 
silenced to zero. For this class of tracer configurations, a 
lengthy but elementary calculation yields 


.F^(Q) — 2{Q X 4- ■ 


(25) 


In particular, if the arms of the cross are chosen to point 
along the directions k = (fc, 0), k^ = (0, k) of the Carte¬ 
sian axes, we have 

X fcj SQ±k±k-^ = -^^Q-i±k,±k) , (26) 

and the above equations simplifies to 


_ ;„2 j-2 J 1, Q = i±k,±k) 




“ 1 0, else 


(27) 


Using this result, we calculate the total energy contribu¬ 
tion of a cross configuration without fast modes as 


E\f,)=2nMk^fl + ^k^ft. 


(28) 


From the MC point of view (24) imposes an extra con¬ 


straint on the allowed phase space in addition to the re¬ 
ality condition /(fc) = f*{—k) for the fast modes during 
the sampling. We thus need to calculate the effect of a 
variation 

Sfiq) = r (dq-fc + -b Sq+k + ) (29) 

of the cross configuration by the real number r on the 
total energy. For the harmonic contribution, it is easy to 
see that 

= ^KNk\rU + rV2) (30) 

It remains to calculate the change of the anharmonic con¬ 


tribution to the energy under a MC move (291. In terms 


of the shift 6T{Q), for which a lengthy and tedious cal¬ 
culation yields 


df{Q) = 2r{Q X kf fiQ - k) + /(Q + k) + r5, 


Q-k-k^ 


+ {k k^). 


(31) 


This last missing piece of information is readily obtained 
from the general variation formula 

2F{Q)5F{-Q) + |<5J'(Q)pj32) 

Q40 

valid for all types of FMC moves. 

Recently [18j we have introduced a new variant of 
FMC that is able to efficiently suppress critical slow¬ 
ing down i.e. exponential growth of integrated autocor¬ 
relation times in critical or nearly critical systems. This 
is achieved by iteratively optimizing the MC acceptance 
rates of individual Fourier amplitudes during the start-up 
phase of the simulation for each wave vector separately, 
aiming at acceptance rates between 30% — 40% for each 
amplitude. In the present simulations, such an optimiza¬ 
tion was, of course, also implemented. 

In what follows, we shall assume that without loss of 
generality k = I, such that only a dependence on the 
anharmonic coupling parameter K remains. The coarse 
graining procedure outlined so far produces a shift 


According to ([^, rescaling of lengths and further 
function” renormalization then leads to 


K I 




k' 

K' 


=: K' 


“wave 

(34) 


since [k] = 0 and [K] = 2. A concrete RG flow K —>■ 
K' is only defined after specifying the function ri{K). 
Imposing invariance k' = I of the harmonic dispersion 
term gives 


r]{K) 


In k{K) 
In b 


(35) 


where we explicitly indicate the dependence of k on the 
parameter K. On the other hand, an invariance condition 
K' = K would implicitly define a function r]i^{K) by 


J ft (Pl ) 

At a FP K' = K = K*, the common value 

r,{K*) = miK*) = rj (37) 


K ■= 






=:K. 


(33) 


of the two functions r]{K) and r]K{K) is the critical ex¬ 
ponent r]. Thus K* can be numerically determined as 
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the location of the common intersection point of these 
functions plotted against K. 

Since the RG transform K i—)■ K'{K) is designed to 
be analytic, we can linearize it around the FP value K* ^ 
such that 


K'{K*+SK)Ki K*+M-6K, (38) 


where 


M := 


dK'{K) 

dK 


K=K* 


(39) 


denotes the slope of the function K K'[K*) — K* 
at its zero K = K*, which can readily be assessed in 
our simulations. For a nontrivial infrared attractive FP 
we expect that K is irrelevant and thus |Af| < 1. If 
K'{K) is not to oscillate back and forth around this FP 
during successive RG iterations, we should also expect 
M > 0. The corresponding Wegner [35] exponent w is 
then defined through AI = b~‘^, i.e. 


InM 


oj = — 


In 6 


(40) 


NUMERICAL EVALUATION STRATEGY 


The membrane systems studied in our simulations may 
be parametrized by a triple of integers (L, I, I'), such that 
A = 2ttI/L,K' = 27rZ'/L,6 = l/l'. Each of the integers 
j = ... ,1' then defines a value 


kj = 2'Kj jL 


for a cross (24) with arms 

kj = (±fcj,0), kj^ = (0, ±fcj) 


(41) 

(42) 


that hosts a tracer configuration with real-valued am¬ 
plitude fc. In principle, each choice of inner cutoff pa¬ 
rameters 1 < I' < I gives rise to fc-values fci,... , fc;/. In 
Ref. [18] we have observed strong finite size irregulari¬ 
ties for the correlation function G{k) = (/(fe)/(—fe)) at 
the smallest accessible nonzero wave vectors, so that it 
is recommended to exclude the fc-value fci from the fol¬ 
lowing fits. Since one needs a minimum of approx. 5 to 6 
values to determine dispersions with sufficient statistical 
reliability (see below), we are confined to lower cutoffs of 
about I' > 6, which in turn puts an approximate lower 
limit of 1/6 > 6/Z on the 6-values accessible in our simu¬ 
lations. On the other hand, due to the discreteness of the 
Brillouin zones of our finite systems, the closest accessible 
value of 1/6 below its ultimate limit 1.0 is 1/6 = {l — l)/l. 

Given such a system, we determine the unnormalized 
probability distribution P^{fc) within a certain interval 
— /max < fc< /max by FMC, and thus the dimensionless 
coarse-grained free energy 

(«) 


To reliably separate the contributions proportional to // 
from those proportional to /^ in a comparison of (28) to 
these data requires to choose a suitable value for /„ 


y max ■ To 

estimate this value, we analyze the bare energy expres¬ 


sion (28), assuming that its coarse-grained counterpart 


will not differ by orders of magnitude from it. After fix¬ 
ing K = 1, our only remaining parameter is the value of 
the remaining bare coupling parameter K. Since both 
the harmonic as well as the anharmonic contribution to 
(28) are proportional to k^, it seems reasonable to choose 
a common value /max for the amplitude at which we ex¬ 
pect to see a factor of A between the bare total energy 
and its purely harmonic part 


for all fc. Numerically, /„ 
tion 


uniformly 

Kn—0 

is determined from the equa- 


Kn 


2KAr/max + ^/max = ^ • 2kn f. 


max 


I.e. 


/max = 2N 


(A — 1)k 


K 


(44) 


(45) 


If A is chosen too small or too large, it becomes numeri¬ 
cally hard to reliably separate harmonic and anharmonic 
contributions by a least squares fit. Moreover, the free 
energy range that needs to be determined in the simu¬ 
lations increases with growing A. After performing var¬ 
ious numerical tests, we settled for a common factor of 
A = 2.6 which was used in all subsequent simulations. 
To actually explore the potential shape in this region, 
a successful simulation approach needs to overcome po¬ 
tentially large free energy differences. After monitoring 
the convergence and tunneling properties of several vari¬ 
ants of the family of multicanonical algorithms, the 1/t 
variant [361137] of the Wang-Landau algorithm [38l [39] 
emerged as a robust and reliable choice. Each single sim¬ 
ulation was performed with an order of magnitude of 200 
tunneling events between /c = 0 and fc = /max at the 1/t 
stage, which took about 10® single MG sweeps per simu¬ 
lation, and for each value of K that we want to inspect, 
V such simulations are needed, one for every fc-value out 
of the set {kj : j = 1,2,..., I'}. An example of the raw 
data obtained by such simulations may be inspected in 

Fig.|3 

The numerical procedure to determine RG recursion 
relations requires therefore a number of nested least 
squares fits. On the one hand, the function E^{fc) ob¬ 
tained from a simulation of the above type will, of course, 
be contaminated by small contributions of powers higher 
than //, i.e. it will not exactly resemble the simple struc¬ 


ture of the bare effective Hamiltonian (13). This well- 


known behavior of generating “new” couplings beyond 
those present in the original bare Hamiltonian, which is 
inherent to the RG, can only be dealt with by fitting to 
a more general ansatz of type 


E'^ifc) = a{k)f/ + 6(fc)// + c(fc)/® + d(fc)/« (46) 















ifc/Nf 


FIG. 1. Raw simulation data for systems of size L = 240 with 
outer cutoff I = 24, inner cutoff I' = 10 and 2d Young modulus 
K — 9.16. The index j = 1, ..., 10 labels the underlying 
tracer configurations of type (241 with stars (421 parametrized 
by the wave numbers kj as defined in Eqn. (41 1 . 


and discarding all coefficients except a{k) and b{k) in 
the analysis that follows. In a second level of fitting, the 
resulting collection of coefficients {a{k),b{k)} is in turn 
fitted to functions of structure 


afit(A:) = 2KNk^ + aefc® + agk^ H-h (47) 

bfit{k) = b^k^ + + • • • + j(48) 


from which the coarse-grained values RntKn are ex¬ 
tracted, while all higher order fit parameters, which cor¬ 
respond to other higher order couplings presumably gen¬ 
erated by the coarse graining operation, are ignored once 
again. 

In practice these fits are not as straightforward to do 
as it may seem. Our numerical tests did show that for 
the first level of fitting, using the truncated eight order 
polynomial (46) as a fit function order gave numerically 
convincing results for all considered parameter ranges. 
However, it turned out to be very hard to decide in ad¬ 
vance how to choose the maximum power Umax of k kept 
in the definition of the fit functions (48) in the second 
level fitting that aims at extracting the lowest order k- 
dependence of the functions a{k) and b{k). Truncating at 


Putting K = 1 without loss of generality, the above pro¬ 
cedure maps every “bare” coupling constants K and shell 
thickness parameter b into a pair of coarse-grained cou¬ 
pling constants {k{K), K{K)). To determine the fixed 
point K* of the underlying RG transformation, it is thus 
necessary to perform simulation for a large number of 


too low orders may yield a certain trade-off among the re¬ 
sulting fit parameters and thus adulterate the results. On 
the other hand, remember that we are necessarily work¬ 
ing in a finite system with a discrete Brillouin zone, the 
minimum spacing between fc-vector components given by 
27r/L, and so only few data points may be available for 
low inner cutoff parameter V. Specifically, for systems 
with small values of, say, I' < 6, a high order polyno¬ 
mial fit will not produce meaningful results, since too 
few unintegrated modes with small fc-vectors parallel to 
a chosen direction are at our disposal. 

Worse, the importance of higher order terms in the 
expansion was observed to strongly vary with the par¬ 
ticular value of K chosen. The “optimal” truncation or¬ 
der Umax of the polynomials may thus even depend on 
K, which makes it very difficult to evaluate the large 
mass of data generated in our simulations in this way. 
Worst, it may be difficult to figure out possible “forbid¬ 
den” powers in the sought-after expansion. For instance, 
we have explicitly checked numerically that no “surface 
tension” contribution cx k^ to a{k) is generated from 
0 by the coarse graining operation, which justifies a 
posteriori the use of 0 as our basic model Hamilto¬ 
nian. Theoretically, this absence can be contributed to 
the presence of a Ward identity (see e.g. the cancella¬ 
tion of /c^-contributions in the sum of contributions to 
Eqs. (30a-d) of Ref. [IS]). However, it is beyond the scope 
of the present work to compute all similar constraints on 
higher order ^-dependent expansion coefficients imposed 
by Ward identities. 

To summarize the above observations, we need to fit 
data obtained for the collection of fc-vectors with a func¬ 
tion of which only the lowest expansion power is known 
with certainty. This problem may look hopeless or at 
least somewhat ill-defined at first sight. Not being aware 
of any pre-assembled approach published in the litera¬ 
ture on numerical mathematics, we had to come up with 
our own custom solution. As is explained in more detail 
in the Appendix, the basic philosophy of our approach is 
not to focus on the unknown higher order contributions 
to the fit functions, but rather to determine the extent 
of validity of the lowest order approximation to the un¬ 
derlying data set and reweighting the members of this 
data set accordingly. Despite currently lacking a rigor¬ 
ous mathematical proof, this seems to work quite well in 
practice. 


different AT-values at each accessible 6 -value to extract 
the values for the fixed point coupling K*{b) and the 
exponents 77 ( 6 ) and w( 6 ) for various system sizes. As 
explained above, for each such iF-value this required I' 
separate Wang-Landau type simulations (one for each kj- 
value) to determine the underlying dispersion changes 
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K 


K 


FIG. 2. Illustration of our numerical p roc edure for parameters L — 240, I = 24 and I' — 10. Left column plots: fits of data 
for k{K) and K{K)/K to the ansatz (491. Each data point shown is derived from fits of the functions a{k),b{k) based on 
up to I' different simulation data of the type shown in Fig. Note the excellent compliance of the numerical data with the 
necessary conditions ii{K) = lirriK^o K{K)/K — 1. Right column, upper plot: illustration of the numerical solution 

of the equations ri{K) = t]k{K) as defined in Eqs. ( |35[ ), ( |36[ ). Right column, lower plot: determination of the slope of function 
K'{K) — K* at K = K*. The direction of the RG flow is schematically indicated. In this particular example, the intersection 
point determined by Eqn. (37 1 is located at K* = 9.80127, and we derive exponents rj = 0.86507 and lo = 1.54456, respectively. 


that govern the computation of exponent rj. For a nu¬ 
merical study it is advantageous to switch to description 
that is continuous in iF, at the same time smoothing the 
statistical noise contained in the resulting data. To an¬ 
alyze the possible iF-dependence of k{K) and K{K), we 
make use of the simple fact that for K —)■ 0 evidently 
both k{K) —)■ k = 1 and K{K)/K —)■ 1, such that any 
analytic fit of these functions must start out with value 
unity at K = Q. Otherwise, the only obvious require¬ 
ment one may impose on a candidate fit function is that 
it should give a smooth and regular interpolation of the 
data. Nnmerical tests have shown that an ansatz of type 

f{K) = 1 -b oln(l -b b'^K) +cK + dK^ + eK^ (49) 


with five free parameters a,... ,e does a good job in this 
respect both for k{K) as well as for K{K)/K (cf. Fig.[^ 
for illustration). With these analytic interpolations at 
hand, it is now an easy task to carry out the analysis 
outlined in Eqs. (35l-(40). 

Error bars for the numerical quantities derived by these 


calculations are derived using a corresponding bootstrap 
analysis [3D] based on 100 bootstrap samples drawn from 
the underlying set of considered coupling values K for 
each particular choice of (L, I, I'). 


RESULTS 

In the present work, we studied three systems with 
one common value of k = 1 and A = 7r/5, with sizes de¬ 
fined by the parameters {L,l) = (120,12), (240,24) and 
(360,36), respectively. Unfortunately, a complete scan 
through all available 6-values was only possible for the 
smallest of these system due to the sheer amount of re¬ 
quired computer resources. The resulting fixed point cou¬ 
pling value K*{b) for L = 120, which is shown in Fig. 
indicates a monotonous fall throughout the whole acces¬ 
sible range. For L = 240, the observed behavior is not in 
conflict with this hypothesis, and we would be very much 
surprised if the behavior for L = 360 were fundamentally 
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lib 


FIG. 3. Numerical results obtained for systems of sizes L = 120, 240, 360 with outer cutoff parameters I — 12, 24, 36, respectively. 
Upper plot: 6-dependent location K*{b) of fixed point value of parameter K. Middle plot: 6-dependent value r]{b) of exponent 
rj. Lower plot: 6-dependent value a;(6) of exponent cu. 


different. Thus, it seems reasonable to assume that there 
is no critical maximum or minimum of K* (b) throughout 
the range 0 < 1 / 6 * < 1 , which only leaves the possibili¬ 
ties 6 * = 1 or 6 * = oo. Give a system defined by (L, I, I'), 
the one with ^ — 1 is of course one for which 1/6 is 

closest to 1. For our present purposes, thin momentum 
shells have some attractive features. For a thin shell the 
number of modes that need to be integrated out during 
the CG step is rather small, such that the simulations 
require less cpu time than for thicker shells. At the same 
time, the large number of remaining unintegrated modes 
inside the shell should increase the numerical reliability 
of determining the dispersions a{k),b{k). Unfortunately, 
however, there is a price to pay for this convenience. In 
fact, for a thin shell, all values k{K) and K{K)/K were 


found to be extremely close to 1 over the whole range of 
considered values of AT, which represents a serious chal¬ 
lenge to a numerical evaluation. For the larger two sys¬ 
tems, the closest accessible estimates for 77 in this limit 
are around r] « 0.93, which is definitely out of range in 
comparison to all other published estimates. Although 
there may exist a common downward trend of 77 ( 6 ) for 
6 —)■ 1 , it is difficult to estimate the limiting behavior. 

Turning to the opposite limit 1/6 —?► 0, we observe 
a nice linear decrease of 77 ( 6 ) with falling 1 / 6 , with all 
the data from various system sizes roughly collapsing on 
a common same master curve, which indicates that for 
determining 77 ( 6 ) finite size effects are small to negligible 
in this limit. Again, it is delicate to extrapolate the data 
to 1/6 —)■ 0. If we assumed that the linear trend persists 
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until 1/6 0 , we would arrive at a rough estimate of 

0.822, (50) 

which is quite satisfying, as it puts our present calcula¬ 
tions roughly in the same ballpark as those done ana¬ 
lytically in the framework of the functional RG 
where the estimate 77 = 0.85 was derived. Nevertheless, 
in view of the finite system size used, the imponderables 
of the above extrapolation 1/6 —>■ 0 and the fact that 
even in this limit the residual error due to the influence 
of irrelevant couplings may only be minimized but not 


completely eliminated, one clearly should not expect (50) 


to be equally precise as the estimate ij = 0.795(10) of 
Ref. [18] derived from a systematic FSS analysis of the 
membrane’s mean squared displacement ((A/)^). 

Extrapolation of uj{b) to 6 —>■ 00 is also delicate. In fact, 
to an unprejudiced reader the data depicted in Fig.|^will 
be compatible with at least two scenarios: 


• Linear extrapolation of oj(b) to 6 —00 produces a 
value of roughly uj = limf,_>.oo w( 6 ) « 4/3 ± 0.3. 

• The 6 -dependence of w( 6 ) may just as well already 
have saturated for 6 —>■ 00 at an asymptotically con¬ 
stant value, leading to any equally crude estimate 
of w « 3/2 ±0.3. 

Even though these estimates may not seem to be ex¬ 
tremely precise, they pave the way for a considerable 
further refinement of our previous FSS result as we show 
next. To explain this in due detail, let us briefly recapit¬ 
ulate the approach followed in Ref. [18] . 

The mean squared displacement ((A/)^) is expected 
to exhibit a finite size scaling behavior of type 

((Aff}r.S + aL^-’>-(l + C(L)) . (51) 


The main obstacle to overcome in an attempt to de¬ 
termine the exponent 77 with high precision is to assess 
the factor ((L) which hosts the subleading corrections to 
scaling by constructing an appropriate ansatz. In princi¬ 
ple, these corrections arise from the presence of irrelevant 
couplings, and thus the leading contributions to ((L) 
should correspond to powers of ..., where 

UJ 2 denotes the Wegner exponent of the next-to-leading 
irrelevant coupling I1III121 (in Ref. [TS] , an additional log¬ 
arithmic contribution of type C,{L) = /31nL±7/L“±... 
has already been ruled out). Unfortunately, however, we 
had (and have) been unable to spot any published nu¬ 
merical estimate for the correction to scaling exponent 
w in the literature. Interestingly, while such an esti¬ 
mate may have been beyond reach for previous simula¬ 
tion approaches to the flat phase of tethered membranes, 
it seems as if w is equally hard to extract from analyt¬ 
ical methods [28]. Lacking any estimate of the correc¬ 
tion to scaling exponent w, in Ref. m we had chosen to 
monitor the error bars for the remaining fit parameters 
produced by different choices of w in the interval [ 0 , 1 ]. 


Based on this reasoning, the estimate 77 = 0.795(10) of 
Ref. [18] quoted above had finally been derived for the 
“naive” choice C{L) = j3/L + 7 /L^. Unfortunately, how¬ 
ever, the sign of the resulting value 6 produced in this fit 
is positive, a fact that was not paid much attention to in 
Ref. [IH]. Recall that according to Eqn. (5) of Ref. [TB] 


((A/)^) ~ J 

where asymptotically for |< 7 | 


'■ <fq 
-> 0 


G{q) 


G{q) = (|/(g)P) 


Kq 


A—rj 


(52) 


(53) 


Based on this ideal power law, a spherical cutoff geometry 
2'k/L <\q\ < K results in 


(A/)^ 


dkl{2TTf _ ± (£/27r)^-^ 


(54) 


I 27 ./L ■' (27r)2(2-?7)K 

which implies a negative value S = — A^“’'/(27r)^(2 — 
r])K, and it is reasonable to expect that this heuristic 
observation also carries over to the case of cubic cutoff 
geometry with subleading scaling corrections included. 
Our present RG approach, in which a;( 6 ) is derived in 


(40) as a by-product of locating the FP coupling K*{b) 


and the exponent 77 ( 6 ) without extra effort, now puts us 
in a position to shed some new light on these problems, 
even though some residual speculations on the structure 
of C(A) are still involved. We propose the ansatz [Kl HT1- 


C(L) = /3/L‘" ± 7 /^"" ± . 


(55) 


which amounts to assuming that either uj 2 ~ 2uj or 
072 3> 2u}, and also to completely discarding additional 
“analytic” corrections of type 1/L. While the first as¬ 
sumption is admittedly difficult to justify based on our 
present knowledge, the latter appears to be reasonable for 
periodic boundary conditions where the renormalization 
of the scaling field gL = 1/L is trivial [45]. To evaluate 
the impact of these scenarios on the numerical estimate 
of 77 , we have carried out new fits of the FSS ansatz ( |^ 
for various choices 0.6 <w< 2.0 to the data for {{AfY) 
generated in Ref. m- The results, which are shown in 
Fig.i reveal a number of interesting points. According 
to Fig. 1^ S < 0 appears to hold only for oj somewhat 
larger than 2 — 77 « 1 . 22 , a range of values for uj that had 
unfortunately not been considered in Ref. [TH]. Our fits 
obviously become singular near this value, but this should 
not come as a surprise: quite trivially, for lu equal to 2 —77 
multiplication of the correction in (55) with 

produces yet another constant besides 5, which results in 
an ill-defined fitting prescription in the close vicinity of 
this value of uj. As the top left plot of Fig. indicates, 
for both our two possible extrapolations w = 4/3 and 
OJ = 3/2, the parameter 6 is indeed negative. We obtain 


77 = 


0.7935(34), w = 4/3 
0.7927(29), w = 3/2 


(56) 
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FIG. 4. Results of least-squares fits of the ansatz ( |51[ ) to 
the data for {{Af)'^) obtained in Ref. |18| for various choices 
0.6 < to < 2.0. Data within the range 4/3 < to < 3/2 are 
plotted in green color as a guide to the eye. Top left panel: 
fit results for 77 . Top right panel: oj-dependence of goodness- 
of-fit parameter Q and reduced x^-parameter of the fits. Bot¬ 
tom right panel: ti;-dependence of fit parameter 5. Note the 
discontinuous change of sign around u) ~ 1.22. Remaining 
panels: cu-dependence of other fit parameters a, P and 7 . 


The accompanying goodness-of-fit parameter Q |4fij dis¬ 
plays a minimum near a; = 3/2 (cf. the top right plot of 
Fig. [^, thus slightly favoring the second of the two sce¬ 
narios (56). In retrospective, we note that their common 
denominator, namely the assertion w > I, could have 
been already anticipated from a FSS analysis similar to 
the carried out in Ref. [18] if only we had monitored the 
sign of the resulting fit parameter 5 in the ansatz (51) for 

m?)- 


In closing this section we note that our new asymp¬ 
totic fitting procedure, which we have used above to com¬ 
pute coarse-grained parameters and which is explained 
in the Appendix, offers yet a complementary way to es¬ 
timate 77 from the FSS of ((A/)^). After all, it was de¬ 
signed for the very purpose of extracting leading func¬ 
tional dependencies from data with unknown higher or¬ 
der corrections. Based on application of the linear ansatz 
log((A/)^) ~ logQ; + (2 — 77 ) logL to the logarithms of the 
data points, the result 77 = 0.7948 obtained by our simple 
recipe is impressively close to that of our above elaborate 
FSS analysis, with nothing more than the leading scal¬ 


ing behavior as input, even though it may be difficult to 
estimate the corresponding error bar. 


DISCUSSION AND OUTLOOK 


In this paper we have illustrated the practical fea¬ 
sibility and usefulness of our FMC implementation of 
Wilson’s MSRG for the flat phase of crystalline mem¬ 
brane, a nontrivial model of continuing physical inter¬ 
est. In particular, we have demonstrated the ability of 
our method to derive the - albeit crude - numerical es¬ 
timate w ± 0.3 S [4/3, 3/2] for the correction to scaling 
exponent in a situation where all other approaches have 
failed so far. Yet, for a meaningful numerical analysis it 
is mandatory to monitor the dependence of observables 
on the thickness parameter h of the employed momentum 
shell. 

Our RG result for cj also allowed to construct an im¬ 
proved FSS procedure for ((A/)^). The resulting new 
estimates (56 1 for the exponent rj deviate even more from 


the value 77 = 0.85 derived both from functional RG 
than our previous one 77 = 0.795(10) given in 
Ref. [18]. In fact, they happen to be much closer to the 
result 77 = 0.78 22(5) extracted from a second order self- 
consistent screening approximation [5D]. On the other 


hand, even our own RG results (50) for 77 show a sim¬ 


ilar tendency to exceed the FSS estimates. A heuristic 
explanation for this common tendency of RG approaches 
to overestimate 77 is as follows. 

For the flat membrane model. Fig. [^indicates that the 
influence of irrelevant couplings is minimal for 6 —>■ 00 . 
Still, even in this limit our RG estimate ( [50] ) is notice¬ 
ably higher than all those obtained from FSS, no matter 
which kind of scaling corrections we employ, even though 
the underlying data were generated using the same un¬ 
derlying FMG algorithm. As Fig. [^ indicates, our RG 
analysis seems not to be afflicted with appreciable finite 
size effects. Thus, the only explanation for this discrep¬ 
ancy is a residual systematic error in the RG result due to 
the influence of the irrelevant couplings that survives the 
limit b —> 00 . In fact, had we not been carefully monitor¬ 
ing the 6 -dependence of our results but had simply cho¬ 
sen one particularly convenient shell configuration, our 
result for 77 might have been still dramatically higher, 
as the middle panel of Fig. indicates. In view of the 
fact that generic functional RG calculations also include 
choosing a projection to a low-dimensional coupling con¬ 
stant space and a cutoff function, these observations may 
hint at the source of the persistent discrepancy between 
RG and FSS estimates for 77 . 

In the near future we plan to investigate the critical 
behavior of hexatic membranes in the so-called crinkled 
phase [TTHSO] using a similar strategy. Gompared to that 
of crystalline membranes, this problem is closely related 
but technically much more involved due to the fact that 
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for hexatic membranes a surface tension contribution of 
type 

J (57) 

has to be taken into account, its coupling constant /i 
being relevant in the RG sense [50]. Using conventional 
FMC we have found it very difficult to determine the 
critical value p,* {k) to which ^ must be tuned at a given 
value of K to actually observe the crinkled phase m- 
Our present approach offers a way to do this, but at the 
expense of determining a two-dimensional flow pattern 
in the variables k and fj, at fixed parameter K. Work in 
this direction is currently in progress. 


Appendix 


To demonstrate the fitting strategy used in extracting 
the lowest order dispersion coefficients from the functions 
a{k) and b{k), we consider a simple toy model defined by 
the function 


fix) = 0.09a;2 - O.lx"^ - 0.2x® -k 0.9x® (A.58) 

in the interval [0,1]. We choose the 20 equidistant points 
Xn '■= 0.05 • n, n = 1,... 20 from this interval, and gener¬ 
ate the telescopic data sets 

Fk ■■= {{Xn, fixn)) : n = 1,..., fc}, fc = 1,..., 20 . 

(A.59) 

We will use the largest considered data set F 20 holding 
20 function values as the input data to our procedure, 
whose goal it is to reconstruct the leading coefficient 0.09 
of fix), based only on the information that fix) should 
start out with a term cx multiplied by a non-negative 
coefficient. 

We start by choosing the fit function 

ifix) := a^x^ (A.60) 


which deliberately ignores all higher order corrections to 
this leading x-dependence that should gradually kick in 
for growing values of x. Of course, the quality of a series 
of least-squares fits of this function applied to the sets Fk 
will successively degrade with growing k. Quantitatively, 
we will observe a crossover from a slow to a steep rise of 
the accompanying y^-parameters 

k 

xl ■= “ ‘P{xi)?, /c = 1 ,... 20 (A.61) 

1=1 


with growing k. The idea is to use the inverse values 

,.2 _ , 

.2 . ^ o (A- 62 ) 


Wk ■■ = 


1/xi, 


Vxfc, 


(since trivially Xi should vanish, we have put wi = 
W 2 in (A.62)) as statistical weights for the data points 



FIG. 5. I llustra tion of the fit procedure outlined in 
Eqs. (A.581-( A.621. Main plot: original data set F 20 pro¬ 
duced from the function (A.581 and final fit with Eqn. (A.601 
(but only shown in the interval (0,0.6)). Inset: normalized 
weights obtained for all 20 data points of the data set E 20 
(note the logarithmic scale). 


(x„, /(x„) in a final least-squares fit of (A.601 to the full 
data set F 20 , which produces our final estimate for a^. 

The described procedure is simple to implement and 
robust, but, of course, far from perfect. Obviously, its 
success relies on a high quality of the underlying data. 
In particular, it is vulnerable to statistical outliers lo¬ 
cated inside the asymptotic region x <C 1. Nevertheless, 
at least in the context of the present paper, where the 
unknown higher order correction terms can be expected 
to be small compared to the leading term, it seems to 
produce sufficiently accurate results. For example, in the 
case of the toy function (A.58), we obtain a = 0.29781, 
which is within 0.73% of the exact value 0.3 = \/0.09. 
In view of the fact that we did not have to make any as¬ 
sumptions on the structure of the fitting function beyond 
its lowest order term, the achieved precision is quite sat¬ 
isfactory. Fig. illustrates the resulting final fit of ipix) 
to the full data set F 20 obtained by reweighting its data 
using the weights Wk successively generated by the pre¬ 
vious fits to the data sets Fk. 
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